Exact results for anomalous transport in one-dimensional Hamiltonian systems 
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Anomalous transport in one-dimensional translation invariant Hamiltonian systems with short range 
interactions, is shown to belong in general to the KPZ universality class. Exact asymptotic forms 
for density-density and current-current time correlation functions and their Fourier transforms are 
given in terms of the Prahofer-Spohn scaling functions, obtained from their exact solution for the 
polynuclear growth model. The exponents of corrections to scaling are found as well, but not so 
the coefficients. Mode coupling theories developed previously are found to be adequate for weakly 
nonlinear chains, but in need of corrections for strongly anharmonic interparticle potentials. A 
simple condition is given under which KPZ behavior does not apply, sound attenuation is only 
logarithmically superdiffusive and heat conduction is more strongly superdiffusive than under KPZ 
behavior. 
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Since the discovery by Alder and Wainwright[l| of 
long-time tails in the Green-Kubo current-current time 
correlations, such as the velocity autocorrelation func- 
tion it has been clear that transport in one and two di- 
mensional Hamiltonian systems must be anomalous in 
most cases. One-dimensional systems have been stud- 
ied extensively in the past decades, both by mode cou- 
pling techniques 0-Q and dynamical scaling and also 
by computer simulations [^H, Most studied are the 
exponents a describing the divergence of the coefficients 
of heat conduction and sound damping with system size 
L as L a , and 8 describing the power law £~( 1-<s ) by which 
the corresponding current-current time correlation func- 
tions decay. For both exponents various values have been 
proposed, with 5 = 1/3 for both heat conduction and 
sound attenuation and a — 1/3 for heat conduction but 
1/2 for sound attenuation being the most common ones 
in recent publications. 

Here I will argue that for generic Hamiltonian systems 
the long time behavior of the dynamics can be obtained 
exactly in terms of the scaling functions obtained by 
Prahofer and SpohnQ for the polynuclear growth model, 
which is in the KPZ universality class. The values of 
S and a mentioned above are confirmed. But also the 
coefficients of size dependent transport coefficients and 
long-time current-current correlation functions are ob- 
tained exactly, as well as the scaling functions describing, 
among other things the asymptotic behaviors of the var- 
ious density-density time correlation functions and their 
Fourier transforms. These results hold in all generality, 
for generic short ranged Id Hamiltonians, from weakly 
anharmonic chains up to mixtures of hard points. They 
establish a rare example of exact results that may be 
obtained for non-integrable Hamiltonian systems out of 
equilibrium. In addition the special conditions under 
which such systems do not belong to the KPZ universality 
class will be formulated simply and sharply, together with 
the consequences for long time and short wave length be- 



havior. 

More specifically, I will discuss classical one- 
dimensional N-particle systems described by a trans- 
lation invariant Hamiltonian with short range interac- 
tions and periodic boundary conditions. Following one 
of the ground-laying papers by Ernst, Hauge and Van 
Leeuwen[9[ I will assume that all slow variables of rel- 
evance for the long time behavior of hydrodynamics 
and related time correlation functions are the long-wave 
length Fourier components of the densities of conserved 
quantities, i.e. particle number, momentum and energy, 
plus products of these. This is a crucial assumption. It is 
not satisfied for most exactly solvable models, which have 
additional slow modes, such as solitonsfioj]. For one- and 
two-dimensional systems the method of EHvL has to be 
generalized somewhat: instead of assuming that the time 
correlation functions of hydrodynamic modes decay ex- 
ponentially with time, one has to write down the mode 
coupling equations as a set of coupled nonlinear equa- 
tions for these correlation functions that must be solved 
self-consistently 0-0] . 

EHvL define the hydrodynamic modes, to leading 
order in the wave number k as linear combinations 
of the Fourier transforms of the microscopic densi- 
ties of particles, momentum and energy (l8|. p M (fc,t) = 
Y,*!=iMf expi-ikxj) - 5 k0 (M(k = 0)), with Mf = 
l,Pj,ej for the particle density n(k,t), the momentum 
density g(k,t) and the energy density e(k,t) respectively 
(l9j . The hydrodynamic modes are two sound modes [20] 
a\{k,t) and a,-i(k,t) and a heat mode an(k,t), given 
respectively, to leading order in k by 

( 8 \ 1/2 

a«r(M)=(j-) (c^KM) + *$(M)) , (1) 



aj?(fc,t) 



nTC z 



1/2 



(e(M)-MM))- (2) 



Here, <r=±l, T is the equilibrium temperature, n the 
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equilibrium number density and p—nm; C p =T(ds/dT) p 
is the specific heat per particle at constant pres- 
sure p, with s the equilibrium entropy per particle; 

1/2 

co=(dp/ dp) s is the adiabatic sound velocity in the limit 
of zero wave number and h is the equilibrium enthalpy 
per particle. Furthermore, 

p(k, t) = (dp/de) n e(k, t) + (dp/dn) e n(k, t), (3) 
7-1 
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-e(k,t) + (dp/dn) e n(k,t), 



where 7 = C p /C v is the specific heat ratio and a = 
— n^ 1 (dn/dT) p the thermal expansion coefficient. The 
allowed values of k are of the form k = ^j^- . To leading 
order in k the hydrodynamic modes are normalized un- 
der the inner product (/, g) = ^ (/*<?), with () a grand 
canonical equilibrium average. 

The time correlation functions of the hydrodynamic 
modes satisfy linear equations involving memory kernels, 



= -icTCokS a (k,t) - k 2 [ dTM a {k,T)S a {k,t-T), 
O t Jo 



words, these functions will assume the forms S a (k,t) = 
exp(— iacokt)T, a (k, t), with T, a (k,t) to a first approxima- 
tion real non-oscillating functions. As a consequence the 
mode coupling contributions to M a are dominated by 
those diagrams in which all vertices are of the type V/ CT 
(but, only in the limit k — diagrams in which the first 
and last vertex are of type V~ a ~ a and all the other ones 
of type VI°~ a also contribute to the leading order). All 
other contributions for at least some time will oscillate 
out of phase with the angular frequency acok of the sound 
mode under consideration. The remaining contributions, 
especially so if described in a coordinate frame comoving 
at the speed of sound, can be identified with the terms 
in a similar mode coupling expansion for the fluctuating 
Burgers equation [l2j], 
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dt 
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Here S a (k,t) — (a a (k,0),a a (k,t)) etc. The mem- 
ory kernels may be expressed through a diagrammatic 
mode coupling expansion as a sum of irreducible skele- 
ton diagrams 11]. These consist of propagators, repre- 
senting stationary density correlation functions Sq(£, t a ), 
and vertices representing the coupling of one propagator 
S(£,t a ) to two propagators S^q,^) and S v (£ — q,t a >>), 
with coupling strength iW^ . For the long time dynam- 
ics only a few of these 27 couplings are important; only 
couplings to two sound modes of the same sign or to 
two heat modes may give rise to long-lived perturba- 
tions, all other combinations of pairs of modes rapidly 
die out through oscillations. From EHvLiQ the relevant 
non-vanishing coupling strengths to leading order in k 
can be obtained as 
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Now a central observation is the following: due to the 
first term on the right-hand side of Eq. (|4]) the sound- 
sound correlation functions will have their weights cen- 
tered around the positions x(t) = x(0) ± cot, in other 



dp(x, t) k dp 2 D d 2 p 
dT~ ~ 2~dx + ~A'dx 2 



dt) 
dx' 



(10) 



with A=Sb(0,0), with the density-density time correla- 
tion function Ss(fc, t) defined, in the limit L — > 00 as 



dxe- lkx S B (x,t)= / dxe~ lkx (p{0 1 0)p{x,t)} 



S B (k,t) 



and Tj[x, t) representing gaussian white noise with 
{r)(x,t)r)(x',t'))=2D6(x - x')6(t - t'). The brackets de- 
note an average over the stationary distribution of the 
density field. This is similar to the hydrodynamic equa- 
tions, but simpler because there is only one conservation 
law. The function Ss{k,t) satisfies an equation similar 
to Eqs. (|4|5p . of the form 



dS B (k,t) 
dt 



dTM B {k, T )S B (k,t-T). (11) 



The mode coupling expansion for this memory kernel has 
exactly the same structure as the set of dominant terms 
for the sound mode memory kernel; all propagators cor- 
respond to the same type of correlation function and all 
vertices carry the same weight factor W, in the case of 
the Burgers equation given by W B = K\A. 

From their exact solution of the polynuclear growth 
modelQ Prahofer and Spohn obtained exact expressions 
for the long time, respectively small frequency behavior 
of the function S B (k,t) and its temporal Fourier trans- 
form S B (k, uj). In the infinite system limit, L — > 00 these 
are of the forml 13] 



S B (k,t) = Afp S (^AkH 2 ) 1 '^ 
S B {k,u) = \/ n .'-x\,.\z fPS 



2K 2 |fc| 3J ^ V(2AK 2 ) 1 /2|fc|3/2 
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with the functions fps and fps defined in Eqs. ((5.3) and 
(5.7) ofRef.0. FromEq.riT] one may obtain expressions 
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for the memory kernel in terms of these scaling functions. 
For the full Fourier transform one obtains 



M B (k,uj) = V2An 2 M, 



ps 



with 



|fc|3/2 



(14) 
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where fp S (w) = J* °° dr exp(iwT)f PS (T 2 / 3 ). 

The corresponding expressions for the long time be- 
havior of the sound modes at non-vanishing k are 

S a (k, t) = exp{-iacokt)fp S ([V2V s t) 2/3 k) , (16) 
M a (k,t) = 2V 2 exp(-iacokt)M PS {k,V2Vst), (17) 

with V„ = W°°. 

Next, I consider the wave number dependent sound 
damping constant r(fc) = 2M a (k, 0) and define the sound 
currents as 
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where Ji{k,t) and Jjj(k,t) are the longitudinal current 
and the heat current Q, denoted by EHvL as J; and J\ 
respectively. Eq. (5.11) of Ref.[8| can now be used to 
obtain the leading small- fc behavior of T(k) and long time 
behavior of ( J a (0, 0) J CT (0, t)) as 




T(k) 
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with T^; Euler's gamma function [2 lj and V s > = W~ a ~ a . 

The leading higher order corrections are obtained by 
replacing in the diagrammatic expansion of the mem- 
ory kernel just one pair of vertices of type V° a by ver- 
tices of type y-"-" or V„ H . One easily shows that 
all these terms add contributions proportional to |fc| -1 / 3 
to T(k) and contributions proportional to t~ 7 / 9 to the 
current-current correlation function. Since there are in- 
finitely many such contributions, there seems to be no 
straightforward way of determining the coefficients ex- 
actly. However, estimates based on the simplest con- 
tributing diagrams can be made 1J] . Further correc- 
tions obtain from terms with 4, 6, • • • vertices of type 
V~ a ~ a or V„ ■ Each of these appears to be of the form 
Cfc _A1 for T(k) and Dt~ v for the current correlation func- 
tion, with C and D constants and /x and v of the form 
A t = 1/3 - J2Z 2 m 3 i 2 /^ and v = 2 / 3 + ET=2 2n i ( 2 / 3 ) 3 



respectively, with m.j and rij natural numbers. Again, for 
each exponent there is an infinity of contributing terms. 

The leading long time behavior of Snik^t) is deter- 
mined in similar way by the sum of all contributions to 
M H (k, t) where the first and last vertex are of type 
and all other vertices are of type V° a , all with the same 
value of cr. These terms do contain an oscillating factor 
exp(— iac^kt), but these oscillations are much slower than 
the oscillations in any of the other terms. Since we have 
to include the contributions to M H of either sign of er, 
we cannot express Sh directly in terms of the Prahofer- 
Spohn scaling functions, but we can do so immediately 
for the memory kernel. A simple analysis yields to lead- 
ing order 

M H (k,t) = 2V£cos(acokt)Mp S {k,V2V s t), (20) 

with Vjj = \ Wff\. For the fc-dependent heat conduction 
coefficient and the heat current time correlation function 
this leads to the expressions 

X{k)=nC p D T {k)=nC p 2 ^^(^-]'\ (21) 
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The heat mode correlation function in Fourier repre- 
sentation is given to leading order by 



S H (k, 



2k 2 V 2 J2 a Mps (k, 



V2V S 



(23) 



From this expression and the asymptotic forms of the 
Prahofer-Spohn scaling functions 8] one immediately 
finds the long-time behavior of the heat-heat correlation 
function for fixed k as 



S H {k,t) = exp [~k 2 D T (k)\t\] 



(24) 



Higher order corrections may be obtained in similar way 
as for the sound modes. 

The dynamic structure factor S(k,u), i.e. the spatio- 
temporal Fourier transform of the density-density time 
correlation function exhibits Brillouin peaks at w = ±cofc 
and a central Rayleigh peak, as usual, but, due to the 
anomalous transport the width and inverse height of the 
Brillouin peaks scale with k as Jfc| 3 ^ 2 @, H[ and those 
of the Rayleigh peak with |fc| 5 / 3 22], in contrast to the 
usual scaling with k 2 . Also, the shape of these peaks is 
not Lorentzian, but is given through the Prahofer-Spohn 
scaling functions as 

s(k,u) = ^R0),M0)) 2 Mfc, w ) 

a 

+ {n(0),a H {0)) 2 S H (k,u), (25) 
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with S a (k,u}) the Fourier transform of S a (k,t). 

All results contained in Eqs. (112II24[) hold in the limit 
L — > oo. For finite periodic systems they apply for 
k = 2im/L and for times or inverse frequencies small 
compared to the sound mode traversal time L/cq. For 
correlation functions at k — 0, like in Eqs. (TT51) and (|2"2")l 
this time range may be extended to a value proportional 
to L 3 / 2 . 

It seems fair to pose that the leading long time dy- 
namics of Id hydrodynamic systems belong to the KPZ 
universality class. The sound-sound correlation functions 
to leading order are identical to the density-density corre- 
lation function of the fluctuating Burgers equation, while 
the heat-heat correlation functions are directly express- 
ible in terms of KPZ memory functions in coordinate 
systems moving at the speed of sound. But notice that 
the correction terms decay only slightly faster with time 
and in most cases will not be negligible up to very large 
times. 

As discussed by Prahofer and Spohn[8] a self consistent 
one-loop mode coupling approximation comes remark- 
ably close to the exact solution to the Burgers equation, 
although there are deviations in the scaling functions of 
up to 10% and not all details of the functional behav- 
ior are captured correctly. Similar results are to be ex- 
pected for the one-loop mode coupling approximation to 
Id hydrodynamics. Some care is required, however with 
using published results. In previous analyses Delfini et 
al.|2j and Wang and Li[4| assumed the sound modes were 
linear combinations of momentum density and displace- 
ment field (equivalent to number density in the absence 
of transversal modes), without contributions from the en- 
ergy density. As can be seen from Eq. ([3]) this is justified 
if (dp/de) n = or equivalently, if C p — C v . This is 
the case for harmonic chains, so it will be a good ap- 
proximation for weakly anharmonic chains. For general 
potentials corrections are needed. 

It has been remarked in several places that the charac- 
teristics of heat conduction and sound dissipation change 
markedly under certain special conditions, such as hav- 
ing an anharmonic nearest neighbor potential symmetric 
in the deviations from the average nearest neighbor dis- 
tance under zero pressure d [3] (this includes the FPU-/3 
model ), o r zero pressure in a system of constrained hard 
points [15J. This can be understood as resulting from a 
vanishing mode coupling amplitude V s . In other words, 
the condition for having non-KPZ behavior quite simply 
and generally is ~ (|^f) s = — 1- It is satisfied indeed for 
the classes of systems mentioned above, but in general 
it does not require any of the criteria quoted above, nor 
the condition C p = C v as conjectured in Rcf.|5]. As dis- 
cussed by Delfini et al. [H , the mode coupling under these 
conditions is dominated by the coupling of a sound mode 
to three sound modes of the same type. Sound damping 
becomes almost normal. In contrast to what is stated 
in Ref.[3| it is still superdiffusive, but only logarithmi- 



cally so (see Ref.[16| for the equivalent case of a growth 
model with leading nonlinearity of cubic order). The k- 
dependent heat conduction coefficient diverges roughly 
as |fc| -1 / 2 and the heat current time correlation func- 
tion decays as £ -1 / 2 , both up to logarithmic corrections. 
The latter must be responsible for the gradual increase 
with time of the exponent 8 for heat conduction from 
roughly 2/5 to 1/2, which has been reported for simula- 
tion results 0, Hj]. 

In stationary states the k a ~ 2 behavior of the Rayleigh 
peak implies a nonlinear temperature profile with, for 
large systems a cusp of form \x — xo| 1_Q (the inverse 
Fourier transform) near a boundary located at xq. Such 
nonlinear profiles have been observed regularly in simu- 
lations, but so far I have nowhere seen mention of this 
simple interpretation. 

Like in the case of the fluctuating Burgers equation, 
the mode coupling equations can also be obtained from 
fluctuating nonlinear hydrodynamic equations [l7|. Since 
their structure remains exactly the same, all the results 
obtained above hold for all one-dimensional systems sat- 
isfying the usual Landau-Lifshitz fluctuating hydrody- 
namic equations with finite transport coefficients. The 
asymptotic long time behavior of density-density and 
current-current time correlation functions is independent 
of the values of the transport coefficients, again like for 
the fluctuating Burgers equation. Obviously, transport 
coefficients in nonlinear transport equations may be fi- 
nite even if the corresponding Green-Kubo integrals are 
divergent. Whether this actually will happen for one and 
two dimensional Hamiltonian systems, as far as I know 
is an open question. 

More detailed derivations of the results presented here 
will be published elsewhere. Besides comparisons to ex- 
isting numerical results, new Molecular Dynamics simu- 
lations will be performed. Applications to quantum sys- 
tems will also be studied. They look feasible but will 
require careful consideration of all quantum aspects. 
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